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I construct predictive models of neutrino mass and mixing that have fewer parameters, both in 
the lepton sector and overall, than the default seesaw model. The predictions are #13 = and 
one massless neutrino, with the models having a Z4, or Z2 symmetry and just one extra degree 
of freedom: one real singlet Higgs field. It has been shown that models with an unbroken family 
symmetry, and with no Higgs fields other than the Standard Model Higgs doublet produce masses 
and mixing matrices that have been ruled out by experiment. Therefore, this article investigates 
the predictions of models with Abelian family symmetries that involve Higgs singlets, doublets and 
triplets, in the hope that they may produce the maximal and minimal mixing angles seen in the best 
fit neutrino mixing matrix. I demonstrate that these models can only produce mixing angles that 
are zero, maximal or unconfined by the symmetry. The maximal mixing angles do not correspond 
to physical mixing, so an Abelian symmetry can, at best, ensure that #13 = 0, while leaving the 
solar and atmospheric mixing angles as free parameters. To generate more features of the best-fit 
mixing matrix a model with a non- Abelian symmetry and a complicated Higgs sector would have 
to be used. 
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The plethora of free parameters is a major reason why the Standard Model (SM) is considered to be only a partial 
ly-j \ theory. Most of the free parameters in the SM are associated with the existence of the three families of fermions. If 
neutrino mass is included, there are six masses, three mixing angles and one phase in the lepton sector. If neutrinos 
are Majorana particles there are an additional two phases in the mixing matrix, and if neutrinos gain mass via the 
seesaw mechanism Q, 0, 0, IH with three right handed neutrinos there are another nine free parameters associated 
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with the heavy right handed neutrinos, giving a total of twenty-one parameters associated with the leptons 00- If 
Oh, only two right handed neutrinos are present the leptons have fourteen parameters. 

A way to reduce the number of free parameters is to add a family symmetry to the SM. A family symmetry is a 
particularly intriguing option for the lepton sector because when neutrinos are massless there is an accidental family 
symmetry J7(1)l c x U(1)i x U(1)l t - Since neutrinos are very light, this family symmetry can be considered to 
be only slightly broken. If there are three massless right handed neutrinos there is also an accidental symmetry 
U(3) UR , which consists of transformations of the right handed neutrinos. It is possible that a subgroup of the 
U(1)l c x U (1) l m x U(1)l t x U (3)j, B family symmetry remains unbroken once the neutrino mass generation mechanism 
is added to the minimal Standard Model. 

A lepton family symmetry has also been suggested as a way to explain the interesting features of the neutrino 
mixing matrix. The mixing matrix, or Maki-Nakagawa-Sakata (MNS) matrix, parameterised by 

C13C12 C13S12 si3e 10 

Uums = I -S23Sl3Cl2e" 10 - C23S12 -S23Sl3«12e~ 10 + C23C12 C13S23 I (1) 
-C23Sl3Cl 2 e" 10 + S23S12 -C 23 Sl3Sl2e"^ - S23C12 C13C23 

(omitting Majorana phases), has two large mixing angles, and one small. The atmospheric mixing angle, #23, is 
maximal at best fit H, 01 El El 0] > an d the small mixing angle, #13, only has an upper bound 0. So, if 613 is 
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taken to be zero, the best fit mixing matrix has the form 
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where 6*12 is the solar mixing angle which is large but not maximal 0, 0, 0, ll^, fl8L HsL l2*oj . This mixing matrix 
pattern may be a result of a family symmetry. If sin #12 = ^7= , Eq. [21 becomes the tri-bimaximal mixing matrix 



|2lL |22j , which hopefully could also arise from a symmetry. 

Many authors have created family symmetry extensions to the SM that attempt to explain aspects of the neutrino 
masses and mixing. Some recent examples are references [2^, l24l I25I I2H . [27j • Common to nearly all the models is 
an extended Higgs sector, involving additional Higgs doublets, triplets or singlets. (The exception is [27j which has 
only the SM Higgs and an explicitly broken symmetry). Many of these models require further complexities, including 
explicit symmetry breaking and relationships between some of the Higgs vacuum expectation values (VEVs), which 
in some models are justified usin g an analysis of the Higgs potential. 

The aim of my previous work |28l I29I ] was to find out whether the mixing matrix of Eq. [21 can be generated by 
a simpler model than the family symmetry models already suggested. The approach I took was to investigate the 
simplest possible family symmetry extensions to the SM, and if these models could not produce the required mixing, 
more complicated models were investigated. 

It was shown in reference [28] that family symmetry models involving just the SM Higgs field can only give neutrino 
masses and mixings that are ruled out by experiment, or give mass matrices that are completely unrestricted. This 
was done by segmenting the study into non-Abelian and Abelian groups. Non-Abelian groups were easily ruled out, 
as they forced at least two of the charged leptons to be degenerate. Abelian groups were investigated by looking at 
all the possible mass matrix textures that can be generated by Abelian groups, and finding their mass eigenvalues 
and mixings angles. 

Since the simplest family symmetry extensions to the SM cannot provide us with a predictive mixing matrix, any 
model that does pred ict the neutrino mixing matrix will require Higgs fields in addition to the SM Higgs doublet. 

The models of j2!l involve a number Higgs doublets that transform under the family symmetry. All possible forms 
of mixing matrices that can be derived from an Abelian symmetry with these Higgs fields were found. #13 = could be 
produced, but the atmospheric mixing and solar mixing angles could not be constrained by the symmetry. Although 
these models reduce the number of parameters in the mixing matrix, the total number of free parameters is larger 
than in models without family symmetries. Also the additional Higgs doublets may allow flavour-changing neutral 
currents. Non-Abelian groups were not studied, as they could not be systematically searched easily. 

This article continues the search to find a model that can generate the neutrino mixing matrix, by extending the 
search to involve the possibility of Higgs singlets which can contribute to a Majorana mass term for the right handed 
neutrinos, and Higgs triplets which can give a left handed neutrino Majorana mass term. Since Dirac neutrinos were 
fully investigated in reference (2i|, this paper assumes neutrinos are Majorana particles. 

This paper is structured as follows. In Sec. [H]the possible Higgs lepton coupling terms are outlined and the family 
symmetries are described. Sec. IIIII shows that for Abelian symmetries, with any number of Higgs fields, the only 
best-fit mixing angle that can be produced is #13 = 0. Sec. IIVI discusses the possibility of using a non-Abelian family 
symmetry to constrain the lepton masses and mixing matrix. This section shows that for models with simple Higgs 
sectors non-Abelian symmetries lead to unrealistic lepton masses. In order to get viable models using non-Abelian 
symmetries either two Higgs doublets are required, or the neutrinos must gain mass using three mass terms: a Dirac 
mass term, and Majorana mass terms for both the left and right handed neutrinos. Sec. searches for the minimal 
model that produces #13 = in the neutrino mixing matrix. Since it has already been shown that models with two 
Higgs doublets can generate #13 = 29], I search for models that have one Higgs doublet, but may have Higgs triplets 
or singlets. Sets of Majorana and Dirac mass matrices that produce #13 = 0, and are compatible with models that 
have only one Higgs doublet are listed. The minimal models that produce #13 — 0, and allow the experimentally 
known values for the other mixing angles, involve the seesaw mechanism and only one real Higgs singlet, in addition 
to the single SM Higgs doublet. The details of these models are presented in Sec. IVII and the number of free variables 
in each of the models is discussed. Sec. I VIII demonstrates that the predictions of these models are unchanged by 
one-loop renormalisation group running. 



II. THE ACTION OF THE FAMILY SYMMETRY 



The neutrino mixing matrix C/mns arises from the diagonalisation of the neutrino mass matrix M„ and the charged 
lepton mass matrix Mf. C/mns = Uj L U u , where Uil diagonalises MiM\ , and U v diagonalises the Majorana matrix 
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M„. The mass matrices arise from the Higgs-lepton couplings, and it is these couplings that are directly affected by 
the family symmetry. Exactly how the leptons couple to the Higgs fields plays an important role in producing the 
mixing matrix via a family symmetry. 

Charged leptons gain mass through a Dirac Yukawa coupling to doublet Higgs fields: 

TXt&itR, (3) 

where Yn are the Yukawa coupling matrices, and L — (vl^l) is the left handed lepton doublet. With only Higgs 
doublets, neutrinos can gain Dirac masses via the term 

LY V (4) 

Majorana masses via a dimension-5 operator |30| from 

~Z$ e iKiJ *fL D , (5) 

and Majorana masses via the seesaw mechanism [E Q IE @ when the right handed Majorana neutrino mass comes 
from a bare mass term: 



(v R ) B M R v R . (6) 

Possible mixing matrices arising from a family symmetry and these three neutrino mass terms have been investigated 
in H3 and Hf. 

Other ways of getting neutrino masses include the left handed leptons coupling to triplet Higgs fields Aj [3l| : 

L\iAiL c , (7) 

where A, are 3x3 coupling matrices. If the neutral components of the triplets gain VEVs, they will generate a 
neutrino mass matrix given by M u = A, (A?). The seesaw mechanism can be altered by adding Higgs singlets Xi to 
give a Majorana mass matrix to the right handed neutrino via 



{v R ) c etXiVR, (8) 

so the right handed Majorana matrix is Mr = e%{xi) + Msare- If the singlets are complex, the neutrinos will also 
couple to the complex conjugate of the singlets x*. Any combination of these mass generation mechanisms can be 
used to give neutrinos mass. The most general light neutrino mass matrix is given by M v — Ml — MoiracM^ 1 Mp irac , 
where Ml is the mass term from the triplet Higgs and/or the dimension-5 operator, Mr is from the bare mass term 
and/or the singlet Higgs coupling term, and Muirac is the Dirac mass matrix due to Higgs doublets. 

The family symmetry alters the mass terms by ensuring each mass term is invariant under the transformations 
which act on the different families of Higgs fields and leptons. The leptons transform via 

L -> X L L, l n -> X m £ R , v R -> X uR ur, (9) 

where the X matrices are sets of 3 x 3 unitary transformation matrices in family space. The families of Higgs doublets, 
triplets and singlets transform via 

$i -> V&jQj, A, -> VAijA, Xi -> V xi jXj, (10) 

where the V matrices are unitary transformation matrices in Higgs family space. The transformation matrices for 
each type of field form a representation of the symmetry group. 

These symmetries dictate the form of the Higgs-lepton coupling matrices, and therefore dictate the form of the mass 
matrices, and hence the mixing matrix pattern. The remainder of this article investigates ways that combinations of 
mass terms and symmetries can generate allowed forms of the neutrino mixing matrix. 



III. ABELIAN FAMILY SYMMETRIES 



To exhaustively search through all the possible family symmetries all groups need to be examined, and all 3 x 3 
matrix representations of these groups need to be considered as candidates for the transformation matrices. However 
there are simplifications that make the systematic study of Abelian groups possible. This section demonstrates that 
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all possible mixing matrix forms that can be created by an Abelian symmetry can be generated from mass matrices 
containing texture zeros, and shows that the only aspect of the best fit neutrino mixing matrix (Eq. [2} that can be 
created by an Abelian symmetry is 6*13 = 0. 

For each group there are many equivalent representations, but fortunately only a finite number of non-equivalent 
representations. Two matrix representations A,; and Yi, are considered to be equivalent if they are related by a 
similarity transformation Xi = S^YiS, where S is a unitary matrix. Appendix lAl demonstrates that two different, but 
equivalent, choices of representation for the lepton transformations yield exactly the same mixing matrix predictions. 
Choosing a different equivalent representation is identical to choosing a different weak basis. Every Abelian represen- 
tation is equivalent to a representation where all the transformation matrices are diagonal, so to study Abelian groups 
only diagonal Xl, X^r and X v r matrices need to be considered. These diagonal lepton transformations either dictate 
that an element of a mass matrix is zero, or the element is completely unrestricted by the symmetry. So the study 
of Abelian groups can be reduced to the study of mass matrices with texture zeros where all the non-zero elements 
are unknown. In fact Grimus et al. have shown that all possible textures of this type can be created from a diagonal 
Abelian symmetry if there are enough Higgs fields in the model [s3 • 

To demonstrate that diagonal representations give texture zeros consider the ijth element of a neutrino mass matrix 
arising from a Higgs triplet. 

M« = A«(A a >, (11) 

where is the zjth element of the Higgs-lepton coupling matrix that couples the ath Higgs triplet to the leptons. 
A symmetry will ensure 

A2' = £*i**W A W J ' (12) 



(no summation over i and j) showing that the matrix element X lJ for a particular Higgs field can only be related to X 13 
for other Higgs fields, and not to other elements of A. Therefore, the ijth element of the mass matrix M v cannot be 
related to any other element (aside from the symmetric relation M V ij = M V ji). So apart from the fact that Majorana 
mass matrices are symmetric, all non-zero elements of M v are completely unrestricted by the symmetry. Similarly, 
Dirac mass matrices for the charged leptons and the neutrinos will get texture zeros from the Abelian symmetry, as 
will the right handed neutrino mass matrix, and neutrino mass matrices due to a dimension-5 operator. 

Note that when the neutrino mass matrix arises from the seesaw mechanism, it is the constituent matrices (-Ml, 
Moirac and Mr) that gain texture zeros from the symmetry. The resultant left handed neutrino mass matrix given by 
M„ = Ml — MuiracM^ 1 Mp irac , may not have texture zeros, and the non-zero elements may be related to each other. 
Most systematic studies of neutrino textures consider the textures zeros in M u only, for example [33l 03 HH |3(| (or in 
the case of |37j the textures of M" 1 ), so Abelian symmetries can recover mass matrix patterns in addition to those 
already studied in these articles. 

To find the form of the mixing matrices that can arise from Abelian symmetries a computerised search was under- 
taken that investigated the mixing matrix predictions of all possible mass matrix textures. A similar search was done 
in [29j, but that search did not allow for textures arising from transforming Higgs triplets and singlets, so it searched 
through a smaller number of mass matrix patterns. 

For a given texture of the neutrino mass matrices and charged lepton mass matrix, sets of matrices were created by 
inserting random numbers in the matrix whenever a non-zero element appeared, and both the neutrino and charged 
lepton mass matrices were diagonaliscd to find the diagonalisation matrices U u and U^l ■ The neutrino mixing matrix 
was then created: ZTmns = Uj L U v . 

This was repeated using the same textures, but different random numbers for the non-zero elements. Then the two 
mixing matrices arising from the same texture were compared. If the two mixing matrices had any identical mixing 
angles then the mixing angle was a result of the symmetry. Many textures can create two or three zero mixing angles, 
so I specifically searched for mixing matrices with 6*13 = and two non-zero mixing angles, and also mixing matrices 
that have mixing angle values that are dictated by the symmetry and are non-zero. This search was performed 
using textures arising from both seesaw and non-seesaw Majorana mass matrices. For the seesaw cases, three right 
handed neutrinos were in general assumed, however all M v matrices arising from fewer right handed neutrinos were 
automatically created. This is because models with two right handed neutrinos give identical M„ matrices to models 
with three right handed neutrinos when the third column of M]ji rac has all zero elements, and Mr is either of an 
upper 2x2 block form, or a diagonal form. Models with one right handed neutrino give identical M v matrices as the 
three neutrino models with two zero columns in M£n rac and a diagonal Mr. 

The results of this search showed that there are four possible types of mixing angle predictions from an Abelian 
symmetry. The mixing angles can be zero, maximal, related to masses, or unrestricted by the symmetry. However, 
the maximal mixing is always unphysical as the two maximally mixed neutrinos are degenerate, as explained in the 
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following paragraph. In models where the mixing angles can be related to the masses, the masses are unrestricted by 
the symmetry, therefore it is not possible for an Abelian symmetry to fully predict any physically meaningful angle 
except for a zero angle. This means that values for the solar and atmospheric mixing angles cannot be derived from 
an Abelian symmetry, though the symmetry may dictate that the angles are related to the mass eigenvalues. 
The maximal mixing angle arises when the Majorana neutrino mass matrix M v has one of two forms. 



M v = 



where the X elements correspond to unknown parameters, is diagonalised by 




(13) 
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(14) 



which indicates the second and third neutrino are maximally mixed. Row and column interchanges of this mass matrix 
also give diagonalisation matrices with one maximal mixing angle. This mass matrix describes one Majorana and one 
Dirac neutrino, or equivalently, one unmixed Majorana neutrino and two maximally mixed Majorana neutrinos with 
mi = — m2. Since Sm 2 — the maximal mixing does not give oscillations. In order to get physical mixing from this 
mass matrix symmetry breaking terms could be added to break the degeneracy, which will also deviate the mixing 
from maximal. 
The mass matrix 




(15) 



has similar properties. It is diagonalised by 



U v = 



cos 9 
V2 



1 

sin ( 

"7T 

cos V 
V2 



le 





cos 9 
sin 9e i 



(16) 



Row and column interchanges of this mass matrix also give a maximal mixing angle. Again the maximally mixed 
neutrinos v 2 and 1/3 have Sm 2 = 0, so maximal oscillation is not possible. 

Even if the degeneracy is broken, the maximally mixed neutrinos cannot correspond to maximal atmospheric mixing 
as no combination of possible U/s and the above U„'s will have the maximal mixing appearing in the correct place in 
the mixing matrix. This is true whether or not 6> 13 = (see |2^| for a more detailed explanation). 

Since neither the maximal atmospheric mixing angle nor a particular allowed solar mixing angle can be generated 
by mass matrices with texture zeros, the only aspect of the best-fit neutrino mixing matrix that can be produced by 
an Abelian symmetry is 6*13 = 0. For a review of other ways to get a small 813 see reference |3gj . 



IV. NON- ABELIAN FAMILY SYMMETRIES 



Since Abelian symmetries cannot generate maximal atmospheric mixing or specify an allowed value for the solar 
mixing, it is worth studying non- Abelian symmetries. This article does not demonstrate which kinds of mixing 
angles can be produced by a non- Abelian symmetry, as I did not find a systematic way to investigate all non- Abelian 
symmetries. However, there are a number of classes of models that will always predict unrealistic masses, and therefore 
cannot produce viable masses and mixings. These classes are described in this section. Future model building searches 
should focus on the classes of models that are not ruled out by their mass spectra. 

If all the leptons transform via an Abelian symmetry, but the Higgs fields transform via a non- Abelian symmetry, 
the model will produce mass matrices that have already been considered in Sec. II I II This is because the lepton 
transformations can always be chosen to be diagonal, and these transformations can only act to force elements of the 
mass matrices to be zero or unrestricted, regardless of how the Higgs fields transform. Therefore these models offer 
no advantages over Abelian groups in predicting the neutrino mixing matrix. 

Reference [2^ considered non- Abelian symmetries with models that have only one Higgs doublet. It was shown 
that if the charged leptons transform under a non- Abelian family symmetry, the symmetry ensures that at least two 
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charged leptons are degenerate, in stark contrast to the known masses which have a large hierarchy. So models with 
only one Higgs doublet are restricted to having the left handed leptons and right handed charged leptons transform 
under an Abelian group. Or equivalently, the sets of transformation matrices {Xl} and {X^r}, must commute. If 
there is more than one Higgs doublet the symmetry may not result in degenerate masses, and may provide a viable 
model. 

This also leaves open the possibility of models with one Higgs doublet that have the right handed neutrinos 
transforming under an non- Abelian symmetry, while the left handed leptons and right handed charged leptons have 
commuting transformations. In this case, if the neutrinos are Dirac particles, the restriction on M„Mj from the left 
handed transformation dictates that the mixing matrix is not allowed by experiment psj , and the right handed neutrino 
transformation will only further restrict the mixing matrix. Models with seesaw neutrinos and no transforming Higgs 
are also ruled out by the left handed transformation alone pgj . 

Seesaw models with one Higgs doublet and transforming Higgs singlets can also give unrealistic results. If we assume 
that there are no triplet Higgs fields and Dimcnsion-5 mass terms are not allowed, then the neutrino mass matrix is 
given by M v = —MoiracM^ 1 M]j irac . The constraints on the Dirac mass matrix are enough to show that three flavour 
oscillation is not possible. Appendix IbI demonstrates that in models where there is only one Higgs doublet, an Abelian 
left handed transformation and a non- Abelian right handed transformation leads to Mnirac having at least two zero 
eigenvalues. The resultant neutrino mass matrix M v will also have at least two zero eigenvalues, corresponding to two 
massless neutrinos, and as dm 2 will be forced to be zero, any mixing between the degenerate neutrinos is unphysical. 

However, if the neutrinos gain mass via the general seesaw mechanism M v = Ml — M Dirac^R ^Dirac ^ ne Ml 
component can provide additional eigenvalues for M u , and may produce acceptable masses and mixing angles. 

Therefore, in order to use a non- Abelian group to generate lepton masses and mixing one of two options is required: 

1. A model with two or more transforming Higgs doublets. 

2. A model where the right handed neutrinos transform according to a non-Abelian group, but all the other leptons 
transform according to an Abelian group, and the neutrinos gain mass via a general seesaw mechanism involving 
Dirac mass terms, Majorana mass terms for both the left and right handed neutrinos and transforming Higgs 
singlets and/or triplets. 

Both of these options require a significant increase in the complexity of the Higgs potential, so it is likely that the 
models may end up having more parameters in total than the SM. 



V. THE SEARCH FOR SIMPLE MODELS WITH 8 13 = 

Section lTTTI showed that Abelian family symmetries can, at best, only ensure that 6*13 = 0, so this section investigates 
the simplest models that make this prediction. Reference [29^ demonstrated that #13 can be produced by models with 
two Higgs doublets and a Z4 symmetry. However, this additional Higgs doublet can cause unwanted flavour changing 
neutral currents, and due to the more complex Higgs potential the models have more free variables than the SM with 
neutrino masses and no family symmetry. To search for simpler models than the two Higgs doublet cases, I looked for 
models that involve only one Higgs doublet, but one or more Higgs singlets or Higgs triplets. Not all textures can be 
created by such a model, in particular the Dirac mass matrices must obey the symmetry restrictions Mi — X L MgXgR 

and Moirac — XjMr>iracX v R. 

There are three ways combinations of diagonalisation matrices with zero or unrestricted angles can give #13 = 
while leaving the other mixing angles non-zero. 

1. U v has an upper left 2x2 block, and Uil has a lower right 2x2 block. 

2. U v is either diagonal, or has an upper left 2x2 block, and Uilzi = 0. 

3. U u \s = 0, and Uil is either diagonal or has a lower right 2x2 block. 

Note that a simultaneous interchange of the rows of Uil and the rows or U v corresponds to an interchange of the 
weak basis leptons, and does not alter the mixing matrix, as U\ L U V remains the same. So stating that U v has a lower 
left 2x2 block and Uil has an upper right 2x2 block is equivalent to stating option one above. 

Using the results of the computerised search, this section finds the mass matrices that can generate a mixing matrix 
with 0i3 = using these three diagonalisation matrix patterns. Symmetry considerations demonstrate that the first 
two options cannot be created using models with only one Higgs doublet, but the third option, with U v having one 
zero, can be generated. 
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A. Two 2x2 diagonalisation matrices 

This diagonalisation pattern requires M„ to have an upper left 2x2 block form, and MiM\ to have a lower right 
2x2 block form. To achieve this Me must have at least one column where both the 2nd row and the 3rd row have 
non-zero elements. To get this from a diagonal transformation restriction of Mg = X^MeXm, X^ must have the form 

X L = diag(e la \e ia2 ,e ia2 ), (17) 

ensuring that rows 2 and 3 have identical patterns of zero and unrestricted elements. This Xj, form also puts 
restrictions on the neutrino mass matrix. 

Looking at the non-seesaw neutrino mass matrices first, the restriction on the neutrino-Higgs triplet coupling 
matrices is /tj = VAijX^KjX^, where the triplets transform as Aj — > VAyAj. If 

(18) 

\C E F J 
then 

e -i( ai+a2 ) B e -2' i a 2£) e ~^2 E j , ( 19 ) 

e —i(a 1 +a 2 )(J e ~2ia 2 p g-2m 2 p 

Note that the phases on B and C are the same, as are the phases on D, E, and F. Therefore, the symmetry either 
forces both B and C to be zero, or allows both B and C to be unconstrained. D, E, and F are similarly linked. So 
M v will be made up of some linear combination of 






(20) 

and could never be of an upper left 2x2 block form. Therefore models with one Higgs doublet and a number of Higgs 
triplets cannot generate an upper left 2x2 block U v and at the same time a lower right 2x2 block for Uil- 

For models with the seesaw mechanism a similar problem occurs. Rows 2 and 3 of Muirac must have identical 
patterns due to the form of Xj,, giving a resultant left handed mass matrix from M v = —Mjji rac M^ M^ irac which 
also has the same pattern after row 2 and 3 and column 2 and 3 interchanges. So the mass matrix also must have the 
form of some linear combination of the matrices in Eq. 1201 and cannot produce an upper left 2x2 block matrix. 

So there is no way that an Abelian symmetry with one Higgs doublet can give us 0±3 — via two 2x2 diagonalisation 
matrices. However if there is more than one Higgs doublet in the model, then an Abelian symmetry can generate this 
pattern p9]|.p|. 



B. Uil with one zero element 

The program outlined in Sec. IIIII generated all Dirac mass matrices made of zero and unrestricted elements, 
diagonalised these matrices, and all the diagonalisation matrices with one zero element were selected. The only Dirac 
mass matrices that can produce such a diagonalisation matrix are 

X \ f X \ 

1, O0X, (21) 

X X J \X XX) 

and row and column interchanges of these matrices. Neither of these matrices can be created from a diagonal 
family symmetry restriction with one Higgs doublet. To illustrate this consider the diagonal transformations X L — 
diag(e lai , e la2 , e IQ3 ) , Xi>r = diag(e 1 ^ 1 , e 1 ^ 2 , e 1 ^ 3 ) . The constraint on the charged lepton mass matrix Mi = x\MiXhr 
reduces to M\ 3 = e lty ^ i ~ a ^M 1 ^ . If M^, Mi%i and M^ are non-zero (as is the case for the matrices in Eq. 121(1 this 
means the phases e i(/32 ~ Q3) = e l ^ 3 ~ a2) = e^ 3 "" 3 ) = 1, and /3 3 = a 3 = (3 2 = a 2 - Consequently, the phase e i(/52 ~ Q2) 
will also equal 1, and there will be no restrictions on the element M122 either, indicating that this element is unknown, 
and is not constrained to be zero by the symmetry Therefore, the above mass matrices that give Uil with one zero 
element cannot arise from a model with one Higgs doublet. 
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C. XI v with one zero element from a triplet Higgs coupling 



The last possibility for finding a mixing matrix with #13 = is from models that have U v with one zero element, 
and UgL either a diagonal or of a lower right 2x2 block form. I separate the study of these models into those that can 
be generated by neutrinos with couplings to triplet Higgs bosons, and those that arise from the seesaw mechanism. 

When neutrinos gain mass from a triplet Higgs coupling, the neutrino mass matrix will be a symmetric matrix with 
elements either zero or unrestricted. The computer-aided search demonstrated that only two of these mass matrices 
textures give a diagonalisation matrix with E/,,13 = 0: 



Both of these matrices can be created by symmetries that also generate MgMj in a lower right 2x2 block, so = 
will result in #13 = when the mixing matrix is constructed. The matrix in Eq. 1221 can be generated by models with 
just one SM Higgs doublet through a dimension-5 mass term or by using the seesaw mechanism with a bare mass 
term. Eq. [221 can also be created using one non-transforming Higgs triplet [2^]. This matrix can be generated by the 
U(l) symmetry L e — — L T , but if the symmetry is not broken the matrix predicts maximal solar mixing, but no 
oscillation since Sm 2 ol = 0, and also predicts v-i to be massless. This symmetry has been considered by many authors 
as a broken symmetry (for examples see |4(j, IU, |43, HE HH HM UM ) ■ However, this symmetry acts only on the leptons, 
and not on any Higgs fields, so the SM Higgs VEV will not spontaneously break the symmetry. The symmetry has to 
be either broken explicitly, or part of a larger extension to the SM involving other fields to do the symmetry breaking. 

The matrix in Eq. 1231 also gives t/,,13 = 0, and predicts a massless third neutrino. In this mass matrix 5m\ 2 and 
the solar mixing angle are related in such a way that the experimental values for these parameters cannot both be 
satisfied. One of the three free variables in the mass matrix supplies a contribution to the atmospheric mixing angle, 
and the other two free variables dictate two masses and the solar mixing angle. The third neutrino is massless. If the 
mass matrix elements are chosen to satisfy the measured dm 2 values, the solar mixing will be very close to maximal, 
and not allowed by experiment. See reference |29| for further explanation. 

This matrix can be created with two Higgs triplets using the symmetry U(1)l c x U{1)l^+l t , where one of the 
triplets has a leptonic charge of L e = 2 and L M + L T = and the other a charge of L e = 1, Lp + L T = 1. 

Hence, family symmetry models with triplet Higgs fields can generate 6*13 = via neutrino diagonalisation matrices 
with one zero element. However, these models have solar mixing angles and mass-squared-differences that are not 
allowed by experiment. 



The seesaw mechanism generates more possibilities for the neutrino mixing matrix form than triplet Higgs couplings. 
As well as neutrino mass matrices that contain zero and unrestricted elements, seesaw mass matrices can have related 
non-zero elements. So the seesaw mechanism with singlet Higgs fields can create mass matrices with no zero elements 
that are diagonalised by matrices with one zero element. Unlike the mass matrices for Eqs. |521and[531 these matrices 
do not overconstrain the solar oscillation. 

There are a large number of combinations of Dirac and right handed Majorana neutrino mass matrices that can 
produce a U v with one zero element. Instead of listing them all, I use some symmetry considerations to omit sets of 
matrices that cannot be produced by a symmetry in models with only one Higgs doublet. 

First of all, if the one zero element in the diagonalisation matrix is to correspond to {7 e 3, MgM\ has to be either 
diagonal, or of a lower right 2x2 block. For Xl = diag(e lai , e 1 " 2 , e zaa ), this means that there must be some Xl 
where a\ ^ 02, otherwise MgM[ would not satifsy this requirement. This restriction eliminates all neutrino Dirac 
mass matrices with the same pattern in rows 1 and 2. 

There are two remaining Dirac mass matrices, that combined with 3x3 right handed Majorana mass matrices can 
give 0i3 = 0: 




(22) 



(23) 



D. U v with one zero from a seesaw mechanism and singlet Higgs couplings 




(24) 
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TABLE I: Right handed neutrino mass matrices that give #13 = when combined with the Dirac mass matrix in Eq. 1241 




M Dirac 2 =0X0. (25) 



Note that M^irac of Eq. [^] has a column of zeros. This means the third right handed neutrino does not couple to 
the left handed neutrinos. In fact the third right handed neutrino is not required to create the desired mixing matrix 
|27l l39j | . So for models with only two right handed neutrinos the Dirac mass matrix 




M Dirac 3 = X , (26) 

in combination with some 2x2 right handed neutrino mass matrices can give diagonalisation matrices with one zero 
element. 

The first Dirac mass matrix (Eq. 12411 can be created from the transformations 

X L =diag(e lai ,e ta2 ,e la2 ), X vR = diag(e lai , e la2 , e"* 3 ). (27) 

This transformation also ensures that MiMj has a lower right 2x2 block form. There are 20 different right handed 
neutrino mass matrices that when combined with this Dirac mass matrix via M v = — M^irac iMr 1 Mliiraci §i ye 
013 = 0. These matrices are listed in table [I] 

The Dirac mass matrix of Eq. 1251 can be generated by the transformations 

X L =diag(e lai ,e la2 ,e ta2 ), X vR = diag(e lai , e M2 , e iai ), (28) 

which again, also ensures that MeMj has a lower right 2x2 block form. Note that the first and the third right handed 
neutrinos transform in the same way, indicating that the right handed neutrino mass matrices must have the same 
form after an interchange of the first and third rows and columns. There are three Majorana matrices that have this 
property and give #13 = when coupled to Eq. 12*51 

/ X X 

X0X\,\XXX\,\X0X\. (29) 




Mr = 

\X X 

The Dirac mass matrix of Eq. and a lower right 2x2 block MgMj can be created by the transformations 

X L =diag(e lai ,e la2 ,e ta2 ), X vR = diag(e lai , e la2 ). (30) 
This mass matrix in conjunction with one of the right handed neutrino mass matrices 

(l*)>«*(ll)> (3D 
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Me 



Symmetry 
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Singlet 
transformation 
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Z 2 





X 
X 



X 






Z 2 




1 
-1 



TABLE II: Models that give #13 = for seesaw neutrinos with couplings to one real singlet Higgs. 



generates a mixing matrix with 6* 13 = 0. 

There are also sets of matrices that produce #13 = from the general seesaw mechanism given by M v = Ml — 
MniracMft 1 Mp irac . However, in all these cases the seesaw part of this matrix (— Mrj>j roc M^ M^ irac ) also gives 
813 = 0, so the addition of the Ml Majorana matrix unnecessarily complicates the model, therefore I have not 
presented these matrices. 

All the matrices listed in this subsection have resultant M v matrices that have no zero elements, but the elements 
of M v are related due to the form of the Dirac and right handed Majorana matrices that they are comprised of. 
Every set of Mr and Mj^irac matrices predict 2/3 to be massless giving an inverse hierarchy of neutrino masses. Using 
8m\ 2 = 8.3 x 10~ 5 eU 2 and 5m\ 3 — 2.4 x lQ^eV 2 H3, the neutrino masses are m 3 = OeV, m 2 = 0.049 eV, and 
mi = 0.058 eV or 0.040 eV. 

This section showed that it is possible for models with one Higgs doublet to create a neutrino mass matrix with 
6 13 = 0. These models all have a neutrino diagonalisation matrix U v with one zero element, and a left handed lepton 
diagonalisation matrix Uil which has a lower right 2x2 block form. Neutrino-triplet Higgs couplings are able to 
generate mass matrices that give #13 = 0, but do not allow the known solar mixing angle and 5m 2 . Models utilising 
the seesaw mechanism and Higgs singlets can create #13 = and give no restrictions on solar or atmospheric mixing. 
These models can have two or three right handed neutrinos. 



VI. THE MODELS THAT GIVE X3 = WITH ONLY ONE REAL SCALAR SINGLET 

None of the mass matrices listed in Sec. IVDI can be created with models with only one doublet and no other Higgs 
fields. The minimal additional field needed to produce any of these sets of matrices is one real scalar singlet. Table 
ITT1 lists the mass matrices that can be created using one real scalar singlet, and the symmetry and transformations 
used. There is one model for each of the three Dirac mass matrices in Ea.s l24l26l The right handed neutrinos gain 
mass through a bare mass term, as well as coupling to the singlet, giving a right handed Majorana mass matrix of 
M R = M Bare + e( X ). 

The mass matrices Mi, M Dirac and Mr in model 3 are the same as the model suggested by Grimus and Lavoura 
[27j • Their model was generated using an L e — — L T symmetry, which was explicitly broken to generate the diagonal 
elements in the right handed Majorana mass matrix. The model listed here instead generates the extra elements in 
Mr using the VEV of a real Higgs singlet as described in the table, but if the L e — — L T symmetry is used, a 
complex Higgs singlet with an L e — — L T charge of 2 will also create the required mass matrices. 

These three models demonstrate that it is possible for a model with family symmetry to reduce the number of 
free parameters in the lepton masses and mixing matrix. However, if the total number of free parameters is more in 
the model than in the SM with neutrino masses, it is difficult to justify the addition of one new Higgs field, so the 
number of free variables in these models should be calculated, and compared to models with the same mass generation 
mechanism, but with no family symmetry. 

The number of physical parameters in the lepton sector with three right handed neutrinos, no family symmetry, 
and no Higgs singlets is twenty-one, corresponding to nine masses, six independent angles, and six phases |y,l3- These 
variables arise from the two Dirac mass matrices and the right handed Majorana mass matrix. A basis where the 
number of free variables is evident is the weak basis where Mi and Mr are diagonal and real. These two mass matrices 
contribute three parameters each. This choice of basis completely defines the right handed neutrino basis, but the 
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TABLE III: Models of table [n] in the basis where Mi and Msare are real and diagonal. All parameters in the mass matrices 
are real. 



left handed leptons can still change basis via a multiplication of a diagonal unitary matrix with three independent 
phases. These three phases can eliminate three free variables in Moirac, ensuring that three of the nine parameters 
are real. The total number of parameters is twenty-one: three from Me, three from Mr and fifteen from M Dirac- 

The seesaw mechanism with two right handed neutrinos has fourteen physical parameters corresponding to seven 
masses, four independent angles, and three phases. The diagonal Mi has three parameters, the diagonal Mr has two 
parameters and Mui rac which is a 3 x 2 matrix has nine parameters. 

The number of variables in the mass matrices listed in table [H] can also be reduced by choosing the basis where Mg 
and Msare are real and diagonal, and by using any extra freedom in the basis choice to reduce the number of phases 
in M Dirac and e. The mass matrices for the three models in this basis are listed in table Hill 

The first model has three variables each in Mg, e, and Moirac and two variables in Msare, giving a total of eleven 
free parameters associated with the leptons, ten fewer parameters than the seesaw model with three right handed 
neutrinos and no family symmetry. 

The second model has three parameters each in Me and Msare, four parameters in e and five parameters in M Dirac, 
giving fifteen parameters in the lepton mass matrices, six fewer than in the seesaw model with three right handed 
neutrinos. 

The third model has two variables each in e and Msare, and three each in Mg and M Dirac, giving ten variables, 
four fewer than in lepton mass matrices in the seesaw model with two right handed neutrinos. 

To see whether these models have fewer parameters overall, consider the Higgs potential for these models. The 
most general Higgs potential for one Higgs doublet, and one Higgs singlet that transforms via \ — > — \ is 

V = C^ty + C 2 (0V)2 + C3X 2 + ^4 + C5X 2 0t0; (32 ) 

where Ci_5, are the five free variables of the Higgs potential, three more than a Higgs potential with only one SM 
doublet. So the first model has seven fewer parameters than the seesaw model with three right handed neutrinos and 
no family symmetry, the second model has three fewer parameters, and the third model has one fewer parameter than 
the seesaw model with two right handed neutrinos. 

All three models have fewer free variables than models with no family symmetry and the same number of right 
handed neutrinos. 



VII. MASS MATRIX RENORMALISATION 



The predictions of — and m„3 = in the models in Sec. IV Dl hold at the seesaw scale, so it is useful to analyse 
the effect of the renormalisation group running to ascertain whether these predictions hold at a low energy scale. 

Under one-loop radiative corrections the Majorana neutrino mass matrix in the charged lepton mass basis at a high 
scale M vmass will evolve to a low scale according to 

M umass — RMvmassR, (33) 

where R = diag(R e , R,,,, R T ), and R e , R^, and R T are numbers relating to the charged lepton Yukawa couplings and 
the energy scales |48|, H3 . 



12 




The transformation M vmasa = Uj L M v Ug L changes a Majorana mass matrix into the charged lepton mass basis, 
where Un, is the diagonalisation matrix for the charged lepton. In the basis where the transformation matrices are 
diagonal, all the models presented give 



cos a sin a | • (34) 
— sin a cos a 

For the matrices of Eqs. [53] and 1231 where the neutrinos gain mass via a triplet Higgs coupling, the pattern of 
zeros and unknown elements in the neutrino mass matrix is unchanged by moving to the charged lepton mass basis. 
The one-loop renormalisation group running does not make zero elements non-zero, only the non-zero elements are 
altered. However since the non-zero elements are in any case unrestricted by the symmetry the renormalisation group 
running does not change any of the predictions of these two matrices. 

For the cases where the neutrino mass matrix arises from the seesaw mechanism, the change to the charged 
lepton mass basis, and the renormalisation group running, can be expressed as a change in the Dirac mass matrix: 

M Dirac mass = RUj L M Dir ac- 



The Dirac mass matrix in Eq. |^has the form 



/ ax 




M Dirac = a 2 | , (35) 
V a 3 

so the renormalised Dirac mass matrix in the charged lepton mass basis is 

' d x 0\ 

^Dirac mass = d 2 I , (36) 

\ a' 3 ) 

where a 1 — aiR e , a 2 = i?/x(«2 cos a — C13 sin a) and a 3 = R T {a2 sin a + 03 cos a). Again, since the texture of the Dirac 
mass matrix is the same before and after renormalisation group running, and the non-zero entries in M Dirac mass are 
free parameters, the predictions are not changed by the one-loop renormalisation group running. The Dirac mass 
matrix of Eq. 1251 in table ITT1 is also unchanged by the renormalisation group. 

Therefore, all of the models that have one Higgs doublet and predict #13 = still predict this after one-loop 
renormalisation group running, and still predict that V3 is massless. The predictions of mass matrices due to triplet 
Higgs coupling fEa.s l2*2l and l23(l also remain the same, so the solar mixing produced by these matrices remain outside 
experimental bounds. If a non-zero #13 is observed, or if the absolute scale of the neutrino mass is found to be larger 
than 0.06eV the three models of table UTI will be falsifcd. 



VIII. CONCLUSION 



Family symmetries may provide a way of explaining the best-fit form of the neutrino mixing matrix, however, it 
was shown that Abelian family symmetries can only explain one aspect of the mixing matrix form: the smallness of 
613 . The minimal model that predicts #13 = 0, and has allowed values for the other mixing angles involves one real 
singlet Higgs, a Z2 symmetry and two right handed neutrinos, although a third right handed neutrino can also be 
accommodated. This model predicts an inverse hierarchy, and #13 remains zero under one-loop renormalisation group 
running. The free variables in these models were counted in order to ascertain whether this model has advantages 
over the SM with neutrino masses. The SM with two right handed neutrinos has sixteen parameters in the lepton 
and Higgs sectors, and with three right handed neutrinos there are twenty-three parameters. The Z2 model with two 
right handed neutrinos and a Higgs singlet has one fewer free variable than the case with no family symmetry, and 
the simplest model with three right handed neutrinos has seven fewer parameters than its equivalent model with no 
family symmetry. So these models not only ensure #13 = 0, they also reduce the number of free variables overall. 

To explain aspects of the mixing matrix other than the smallness of #13, a non- Abelian symmetry will be required, 
however in many simple models a non- Abelian symmetry guarantees unrealistic lepton masses. To produce a viable 
non- Abelian symmetry model either more than one Higgs doublet is required, or the neutrinos must gain mass through 
two different mechanisms: The seesaw mechanism and coupling to a triplet Higgs or via a dimension-5 operator. So 
either we have to be content with family symmetry models predicting only a small aspect of the mixing matrix, or 
we have to accept the consequences of more complex models. 
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APPENDIX A: EQUIVALENT REPRESENTATIONS OF THE LEPTON TRANSFORMATIONS 

This appendix shows that two different, but equivalent, representations for the lepton transformation matrices give 
exactly the same restrictions on the lepton mixing matrix. It follows the structure of Appendix A in reference pflj but 
generalises that result. The method is to find the restrictions on the mass matrices due to the transformations, and 
since the transformation matrices are related via a similarity transformation, the restrictions on the mass matrices 
from both sets of transformations are also related. The relationship between the diagonalisation matrices from the 
two representations are then found, and the mixing matrices are compared. 

Two equivalent representations for the fermion transformations, Al, Am, A v r and Bl, BgR, B^r, are related by 

B L = SlA L S L , (Al) 
Ben = S iR AenSiR, (A2) 
B v r — S vR A v rS v r. (A3) 

where Sl,S^r and S v r can be any 3x3 unitary matrices. The Higgs doublets, singlets and triplets transform 
respectively as 

$ -> Vi$, x -> V X S, A -> Va A. (A4) 

Reference j2£j demonstrated that the charged lepton mass matrix derived from the A transformations, Mia, is 
related to the charged lepton mass matrix from the B transformation, Mg, B . Mia has the same restrictions from the 
symmetry as SlMi B S\ r , so the two matrices have equivalent forms. Diagonalising these two mass matrices shows 
that Uf L has the same restrictions from the symmetry as S\uf L , where U^£ B are the left diagonalisation matrices for 
Mia and Mib respectively. Similarly, the Dirac neutrino mass matrix Moirac A has the same restrictions due to the 
symmetry as ShMmmr. bS.. h . 

Reference |29j also shows that neutrino mass matrices arising from Dimension-5 operators or the seesaw mechanism 
with no Higgs singlets have similar properties. M vA has the same restrictions as SlM„bS~[, so by setting M u a = 
SlMvbSJ^ and diagonalising, the diagonalisation matrices can be shown to be related by U„ — S\U^. Putting the 
neutrino and charged lepton diagonalisation matrices together demonstrates that the mixing matrices are the same: 



Next consider neutrinos that gain mass by coupling to a Higgs triplet. In this case, the neutrino mass 
given by M v = Ai(A°). The coupling matrices \ are constrained by the A transformation according to 

^Ai = A\v&ij\AjA* L , 
and are constrained by the B transformation according to 



Rearranging shows that 



tjB _ jjBUjB _ jjAj a of ttA 

- TJ A 

— U MNS- 



\Bi — B\Vj\ij\BjB* L 

= SlAlS L V AtJ \ Bj SlAlS* L . 

{S L X Bi Sl) = A\V Mj {S L \ Bj Sl)Al, 



(A5) 
matrix is 

(A6) 

(A7) 
(A8) 

(A9) 



demonstrating that Sl^BiS^ has the same restrictions due to the symmetry as Xau so they can be equated. Forming 
the mass matrix for A shows 

M vA = A Al (A°) (A10) 

= S L X Bl (A° t )Sl (All) 

= S L M vB Sl, (A12) 
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demonstrating that M v a has exactly the same restrictions due to the symmetry as SlMvb&l ■ Diagonalising shows 
that the diagonalisation matrices must be related by U,f = S L U^, and using Eq. |H the mixing matrices from the A 
and B transformations can be shown to have identical forms. 

When singlet Higgs fields are present, and the right handed neutrinos gain mass via coupling to the singlet Higgs 
and a bare mass term, the right handed Majorana mass matrix is given by M R — ti{xi) + M Bare . The coupling 
matrices ti are restricted by the B transformations by 

e Bi = B^ R V X ij€BjB V R (A13) 

= Sl R A^ R Sl R V X i j e B jSl R A l , R S vR . (A14) 

Rearranging this equation shows 

( s tR e BiSl R ) = A T vR V x ij(S* R e B jS vR )A„ R , (A15) 

demonstrating that S* R e B iS] /R has the same restrictions as CAi- Similarly, the bare mass matrix from the B trans- 
formation, MsareB, is related to that of the A tr 
handed mass matrix from the A transformation is 



formation, M BareB , is related to that of the A transformation by M Bare A = S* R M Bare bS vR . So the total right 



Mr A = e Ai {Xi)+M Bar eA (A16) 
= S: R (e Bi ( X i) + M Bare B )Sl R (A17) 
= S; R M RB St R . (A18) 

Since the two Dirac mass matrices are related by MoiracA — ShMuirac bS\ r , the resultant left handed neutrino 
mass matrix is 

M vA = M DlracA M R \M T DtracA (A19) 
= S L M DiracB Sl R S^RM R 1 B S^ R SlM'[ lir . acB S^ R (A20) 
= S L M vB Sl R . (A21) 

Again, the relationship between the mass matrices means that = S\u£~ , and using Eq. IA5I the mixing matrix is 
the same for both the A and B transformations. 

All methods of gaining Majorana neutrino masses show the relationship M v a = SlM vB SJ i , so a model that has 
neutrino mass gained from a number of these methods will also display this relationship between the mass matrices 
from different, but equivalent representations. This means that for models with any combination of the methods 
shown above, two equivalent representations will still give identical results. 



APPENDIX B: NON-ABELIAN RIGHT HANDED NEUTRINO TRANSFORMATIONS 

The following shows that in models with one Higgs doublet, left handed neutrinos transforming via an Abelian 
group, and right handed neutrinos transforming via a non-Abelian group always result in M£> irac having at least two 
zero eigenvalues. 

For the right handed neutrinos to transform via a non-Abelian group there needs to be at least two transformation 
matrices that do not commute, two of which I label X uR \ and X uR 2 ■ The left handed transformations that act in 
conjunction with these right handed transformations are labelled X^i and X^i. The Dirac mass matrix must obey 
the symmetry restrictions 

M£)i rac = X L1 M£ii rac X vR i = X^ L2 Al£>i rac X,y R 2- (Bl) 

If the basis where X^i, X^i and X vR \ are diagonal is chosen, the first transformation acts to enforce some texture 
zeros on the mass matrix. The only mass matrices that can arise from such a transformation, and have at least two 
non-zero eigenvalues are 



^00 0\ /0 0\ /0a 3 0\ 

a 2 , a 3 , a 2 , 
\00ai/ \ 02 oi / \0 flj 



(B2) 
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and row and column interchanges of these matrices. 

Applying the second transformation to these already constrained matrices will always reduce the number of non- 
zero eigenvalues to zero or one. To illustrate this just consider the first matrix in Eq. IB2I To create this matrix the 
first transformation matrices have the form 



X L i = diag{e la \e la \e la3 ), 
X um = diag(e iai ,e ia *,e ias ), 



and the action of the second transformation is 



M. 



Di 



~ ^L2^DiracX u H2 



( e~ il31 



\ 




where U = X V R2- The restrictions on elements 22 and 33 of matrix Moirac are 



Oi 



aiU 2 2, a 2 



-i/3 s 



aiU; 



33 • 



(B3) 
(B4) 



(B5) 
(B6) 

(B7) 



If both these elements are to remain non-zero, e~ % ^Ui2 = e~ l ^ 3 U^ = 1, therefore | C/22 1 = | C/33 1 = 1, and the unitary 
matrix must be diagonal. A diagonal X U R2 will commute with X^m, meaning that the right handed neutrinos actually 
transform via an Abelian group. This is the same for all of the other matrices in Eq. IB 21 - either the mass matrix has 
two or three zero eigenvalues, or the right handed neutrinos transform according to an Abelian group. 
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